Higher order parametric level statistics in disordered systems 
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Higher order parametric level correlations in disordered systems with broken time-reversal sym- 
metry are studied by mapping the problem onto a model of coupled Hermitian random matrices. 
Closed analytical expression is derived for parametric density-density correlation function which 
corresponds to a perturbation of disordered system by a multicomponent flux. 



PACS number(s): 05.45. +b, 71.55.Jv, 73.23.-b 

Parametric level statistics reflects the response of the 
spectrum {E n } of complex chaotic systems to an exter- 
nal perturbation. A few years ago it was shown 
that a system whose spectrum follows closely the univer- 
sal fluctuations predicted by the random matrix theory 
should also exhibit a universal parametric behavior. 
This conclusion was reached by analyzing the dimension- 
less autocorrelator of level velocities of electron in a dis- 
ordered metallic sample with a ring topology, enclosing 
a magnetic flux (p. A diagrammatic perturbation tech- 
nique was used in the range of fluxes g^ 1 ^ 2 <C tp <C 1 |Q, 
while the opposite limit, ip < g -1 / 2 0, has been treated 
within the framework of the supersymmetry formalism 
HIJ. [Here g 3> 1 is the dimensionless conductance]. In 
this particular problem, the parametric correlations take 
a universal form involving the rescaled parameter X 2 = 
Airgtp 2 , with g — E c /A, the ratio of the Thouless en- 
ergy and the mean level spacing. Numerical simulations 
have supported the point that the universal character 
of parametric level statistics extends to a wider class of 
chaotic systems without disorder [chaotic billiards] whose 
Hamiltonian depends on some external parameter x. In 
such systems, the spectral fluctuations taken at different 
values of x become system-independent after rescaling, 
x — > X, which involves solely the 'generalized' dimen- 
sionless conductance g = (47T)" 1 A~ 2 ( (dE n (x) /dx) 2 ). 
Along with the diagrammatic technique and the super- 
symmetry formalism, the parametric correlations have 
been studied in detail within the model of Brownian mo- 
tion 0^] and in the semiclassical limit |9|JlO|] . 

A further burst of activity in the field occurred after it 
was realized [p^Jll|-[p4]| that the problem of parametric 
level correlations is identical to the ground state dynam- 
ics of the integrable many-body quantum model known 
as Calogero-Sutherland-Moser [CSM] system. This gave 
new information about CSM space-time (r, r) correlation 
functions that can be obtained from parametric density- 
density correlators ( v(E, Q)v{E' , <p) ) involving only two 
different external parameters by mapping || X 2 — > —2ir, 
E/A — > r. For the more general situation of higher or- 
der correlation functions the connection between CSM 
fermions and quantum chaotic systems has been estab- 



lished as well by using the supersymmetry technique, 
however it has not led to any explicit analytical results 
beyond the two-point correlators due to enormous in- 
crease of number of entries in the supermatrix fields, 
thereby making any explicit calculations in that approach 
impossible. Extensions to higher order statistics can be 
performed by using an involved method of differential 
equations for quantum correlation functions proposed in 
the much earlier work |l5| ]. 

In the present paper we address the issue of higher or- 
der parametric level statistics within the framework of 
the random matrix theory, by appealingto the model 
of coupled Hermitian random matrices ||16|] . The latter 
enables us to provide a complete information about para- 
metric correlations of single electron level densities in the 
presence of the multicomponent flux perturbing a disor- 
dered system, characterized by a dimensionless conduc- 
tance g ^> 1. To the best of our knowledge, this is the 
first detailed study of higher order parametric level statis- 
tics in disordered systems which adopts the conventional 
language of the random matrix theory. 

In what follows we consider a weakly disordered system 
fallen in the universal (metallic) regime, j » 1, which 
is known J5| to be modelled by invariant ensembles of 
large random matrices. Assuming that the time reversal 
symmetry is completely broken (unitary symmetry) , one 
can statistically describe an unperturbed single electron 
spectrum by Gaussian Unitary Ensemble [GUE] of large 
N x N random matrices Hq distributed in accordance 
with the probability density V [Hq] cx exp { — Tr H 2 } . 
Such a distribution V [Ho] induces the energy scale A 
being the mean level spacing, A = n (2N) 1 . Let us 
now apply a Gaussian perturbation consisting of d com- 
ponents ~0 d — (</>i,..., <f>d) which does not change the 



global unitary symmetry, and which drives the Hamil- 
tonian Ho to H = Hq + y^u—i 4>kHk, with matrices 
Hk drawn from GUE: 



V [H k ] k exp {- Tr H 2 } for k = 
I, . . . , d. This choice corresponds to the equal 'strength' 
of each component of the 'vector' perturbation Tp* d since 
the average {{H^^viHk)^^') is independent of the in- 
dex k. The quantity which provides the most detailed 
information about parametric correlations in the case of 



f 



the multicomponent perturbation lft d is the correlator of this reason, we will concentrate on the dimensionless mul- 
level densities v(E 7 ~0 a ) — Tr 6(E — H — J2k=i ^kHk) tipoint correlator 
taken at both different values of energy E and of a. For 

/ PO Pi Pd \ 

k PB _ Pd (w o) }^;.--;w d) }^ d )=^ m ( n^+^ n o) , ^ n^+^^i) ••■ n K^+^f,^)), a) 

\i = l ii=l id = l I 

where m = p + . . . + p dl and the angular brackets stand for averaging over ensembles of Hermitian matrices Hf. 
with k = 0, . . . , d. Equation (|l]) can be rewritten as a (d + 1) multiple matrix integral over matrices H = H and 

*po,...,w cx I dH Q ... I dH d HY[Tr6(E + J° ] - cxp I - Tr 

(2) 
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with 0o — l a n( i = oo. [This convention is relaxed 
everywhere below Eq. (^)] . We notice that the strengths 
(/)). (k = 1, . . . , d) of the perturbation are supposed to be 
small, <fik *C !• This is justified in the thermodynamic 
limit 2V — > oo, since for Gaussian perturbation accepted 
above, 0fc are known to scale with N as <f>k = nN~ 1 / 2 Xk, 
with Xfe being the set of dimensionless parameters of or- 
der unity Q . 

Our crucial observation is that Eq. (^) can be in- 
terpreted as a density-density correlator in the effective 
model of (d + 1) Hermitian random matrices coupled in 
a chain: Each matrix H a is represented by a point, and 
two adjacent matrices H a and H a+ i are joined by a line 



if the coupling of the type exp{c Q Tr H a H a+ {\ is present 
in Eq. (g). In this situation, the joint probability den- 
sity of eigenvalues of all the matrices in the chain can be 
deduced through the Itzykson-Zuber integral jL8] mak- 
ing the model of random Hermitian matrices coupled in 
a chain to be a completely solvable. In accordance with 



the Eynard-Mehta theorem |16|, the dimensionless cor- 
relator fcp 0j ... ! p d can be represented as a determinant of 
the to x to block matrix, m = pa + . . . + Pd, consisting 
of (d + 1) x (d + 1) rectangular submatrices K a ^ with 
a, f3 — 1, . . . , (d + 1), each of them having p a -\ x pp-i 
entries fn9|, 
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The matrix kernels AT Qj/ 3 in Eq. (||) are 

(£, 77) = A [flo,/s & »?) - E aS (£, 77)] 

where 



iV-l 



(4) 
(5) 

E a ,p (£,v) = ( w a * ■ ■ ■ * wp-i) {£,rj) (6) 

for 1 < a < < d + 1; otherwise, E a< p = 0. Here the 
partial weights w a are 

v a (0 + v a+1 ( v ) 



and 



Wa{£,ri) = exp{- 



+ 20"^}, (7a) 



v a (0 = (C-i + K z )K,i + + (7b) 



(compare with the weight of the matrix model, Eq. (0)). 
The notation (w a * ... * tf/3-i) 77) stands for the prod- 
uct of the partial weights w integrated over internal vari- 
ables of that product. Two sets of orthogonal functions 
P a j and Qp.j entering Eq. (||) are determined recursively 



p ».j (£) = / dr)P a -i,j (77) (77, () , (8a) 
Qpj (0 = / dl l w P (^ 7 ?) Q 13+1,3 (v) , (8b) 



for 2 < a < d+l and 1 < (3 < d; the starting points of the 
recursions (pal) and (pB) are the polynomials Pij = Pj 



2 



and Qd+ij — Qj orthogonal with respect to a nonlocal 
weight W (£, rj) = (wi * . . . * w d ) (£, 77), 



(9) 



Close inspection of the equations above shows that 
the basic orthogonal polynomials Pj and Qj can be ex- 
pressed in terms of Hermite polynomials, Pj(£) = Hj(£), 

Qj(0 = Hj(£[l + ELi^fc] -1/2 )' Then, step-by-step 
integrations in Eqs. (jg) yield 



x exp 

for f3 > a + 2, while E a , a - 
the function G a reads 



and _E QjQ ,+i 



Ga (0 



2 



(14) 
Here 

(15) 



Qa,j (0 



J/2 



C Q -i 



1 + 



, (10a) 



(10b) 



where we have introduced the Hermite functions ]2C| ] 
$j (£) = exp [-£ 2 /2] flj (0- A1 so, we defined the func- 
tion 



(11) 



and the constant C a = [1 + 2fe=i ^fc] 1 ^ 2 - P n order to 
compactify the formulas, it is agreed from now on that 

4> d+l = <p d , 00 = 4>i, Tiz (•••) = o, and nf<: (■ •_) = 

1]. One can verify that the orthogonality relation (0) is 
satisfied with 



fc=l fc=l 
so that the first term in Eq. (0) is 



(12) 



JV-l 
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(13c) 



The second term in Eq. (0) is found from Eqs. and 
(0), 



Now, we are in position to compute the matrix kernels 
K a ^ via Eqs. (||), (jl^) and ( |l4] ) in the leading order in 
N ^ oo and keeping X fc = ^N 1 ' 2 /^ ~ O (1) fixed. The 
simplest, diagonal kernel K a ^ a can be evaluated through 
the Christoffel-Darboux formula | pl[ , supplemented by 
the asymptotics of Hermite functions, 



$2W (t) 
<&2Af+l (*) 



(-1) 



N 



cos {2tN 1 ' 2 ) 
sin {2tN 1 / 2 ) 



where t ~ AO (iV°) . One obtains, 

_ P G„(a-cfa) sin t 7 ^" 1 (£ ~ *?)] 



(£, J?) 



(16) 



(17) 



Two other cases, a < (3 and a > (3, demand more effort. 
For a < (3 we represent the sum for H a< p in Eq. ( |l3b| ) 
as a difference of two series, (• ■ — E^Af (••■)■ 

The first sum is exactly computable by making use of 
the Mehler summation formula pi] ]. In the thermody- 
namic limit, this procedure yields a term which is equal 
to E aj p in Eq. (O ), and therefore it gets canceled from 
the expression ([|) for K a< p which is completely due to 
the remaining sum Y^Ln (■■■)■ To evaluate the latter, 
we replace the sum over j by an integral to get 



K, 



a</3 



X / d\i cos|7r — — 
'i A 



2\2 0- 1 
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-Ai}exp{-- 



X 2 }. (18) 



In the case a > f3 the large-j terms in Eq. (13c) 
yield the main contribution to the sum due to the fac- 
tor [1 + Ylk=d 0/cP^ 2 - Then, passing from summation to 
integration, we derive 



K. 



a>0 



i 



d\ cos{7r — - — A} exp{ 



9i9 a — 1 
?7,, r 7T A 



fc=/3 



(19) 



Notice that the structure of the block matrix in Eq. (|^) 
allows one to simultaneously suppress the prefactors of 



3 



the form ]J k (• • ■) ^ in Eqs. |t7|), © and @. Hav- 
ing this in mind, we come down to the closed analytical 
determinantal expression Eq. (||) for (po + . . . +pd)-point 
density-density correlator with K a ^ replaced by M Qii a, 



M a<f3 (£, i]) 



_ sin [ttA- 1 (g - 77)] 
ttA-i (£ - ry) 



/oo 
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(20a) 



(20b) 



M. 



a>/3 ' 



fan) = I dAcos{7r^A} 



7T 2 X 2 

fc=/3 



(20c) 



Equations (||) and (|20| ) are the main result of the pa- 
per. They provide a detailed information about higher 
order parametric density-density correlations in the case 
of multiparameter perturbation of disordered system. 
Several particular correlators can be readily deduced 
from our general expression: (i) For the scalar per- 
turbation, one obtains that k Ptq = A p+q ( Y[ P i=i V {E + 
) is determined by 



.o)nLi"(^ + nj 



k Pt q = Det 



M ata (fA}i,(jJj) M a< p{Wi,Vlj) 

M a>p (Qi,Uj) M a , a (Oi,%) 



(21) 



where M a ^ are those given by Eqs. (20) with J2k 



A 2 ; (ii) By replacement ^ uJi/A — > r», f2j/A — > i?i and 
X 2 — > — 2ir in Eq. ( pl| ) one arrives at the space-time 
correlation function in the CSM model with a coupling 
A = 1; here the coordinates {r^} correspond to the time 
t = 0, while the {Ri} refer to the time t = r. 

In summary, we presented a random-matrix-theory 
treatment of the problem of higher order parametric spec- 
tral statistics in disordered systems with broken time 
reveral symmetry in the presence of the multiparameter 
perturbation. A complete analytical solution was based 
on the mapping the initial problem onto a model of ran- 
dom Hermitian matrices coupled in a chain. As a partic- 
ular case of the general solution given by Eqs. (||) and 

(20) , the multipoint parametric spectral correlator Eq. 

(21) for the scalar perturbation has been obtained. To- 
gether with a well established correspondence between 
CSM fermions and parametric level statistics, the latter 
expression provides an information about the space-time 
correlation function in the Calogero-Sutherland-Moser 
model of free, non-interacting fermions. 

The authors thank I. Yurkcvich for bringing the refer- 
ence [ll5| to our attention. 
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